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According to the Heisenberg uncertainty principle of quan- 
tum mechanics, particles which are localized in space by a 
bounding potential must have a finite distribution of mo- 
menta. This leads, even in the lowest-possible energy state, to 
vibrations, and thus, to the so-called zero-point energy. For 
chemically bound hydrogen the zero-point energy can be quite 
substantial. For example, for a free H2 molecule it is 0.26 eV, 
a significant value in the realm of chemistry, where often an 
energy of the order of 0.1 eV/atom (or 2.3 kcal/mol) decides 
whether or not a chemical reaction takes place with an appre- 
ciable rate. Yet, in many theoretical studies the dynamics of 
chemical reactions involving hydrogen has been treated clas- 
sically or quasi-classically, assuming that the quantum me- 
chanical nature of H nuclei, i.e. the zero-point effects, will 
not strongly affect the relevant physical or chemical proper- 
ties. In this paper we show that this assumption is not justi- 
fied. We will demonstrate that for very basic and fundamental 
catalytic-reaction steps, namely the dissociative adsorption of 
molecular hydrogen at transition metal surfaces and its time- 
reverse process, the associative desorption, zero-point effects 
can not only quantitatively but even qualitatively affect the 
chemical processes and rates. Our calculations (treating elec- 
trons as well as H nuclei quantum-mechanically) establish the 
importance of additional zero-point effects generated by the 
H2-surface interaction and how energy of the H-H stretch vi- 
bration is transferred into those and vice versa. 



Recently density functional theory calculations of the 
important diversity of pathways for chemical reactions 
at extended metal surfaces became possible 0-^. Still, 
these calculated high-dimensional ab initio potential- 
energy surfaces (PES), which fully determine the energet- 
ics of the nuclei, are too complex that their consequences 
on the chemical reaction dynamics were obvious. It is 
worth to stress at this point that these calculations re- 
vealed that the high-dimensional nature of the PES can 
be very important, i.e., a restriction to just one reaction 
pathway is often not sufficient ^g- In order to analyze 
and understand the catalytic dissociation it is important 
and possibly even necessary to combine the information 
about the high-dimensional PES of the chemical reaction 
with a subsequent calculation of the dynamics of the nu- 
clei, moving along this PES. Typically this dynamics is 
described applying the laws of classical mechanics, the 
main motivation of this approach being its simplicity. 
Some studies of dissociative adsorption (and associative 
desorption) of H2 , which is the topic we concentrate on in 
this paper, had addressed the issue of comparing classical 



and quantum dynamics IQ-^ , but in these studies only up 
to four degrees of freedom of the hydrogen molecule were 
treated quantum mechanically; obviously, the dynamics 
of two nuclei is determined by six coordinates. In the fol- 
lowing we will restrict our discussion to the dynamics of 
dissociative adsorption. However, our findings about the 
importance of zero-point effects are equally valid for the 
associative desorption since the processes of dissociative 
adsorption and associative desorption are related by the 
principle of microscopic reversibility. 

In this paper we will show that a quantum-mechanical, 
six-dimensional treatment of the nuclei exposes several 
interesting and important properties. Such calculations 
became possible only recently Q, and as a first ex- 
ample we studied the dissociative adsorption of II2 at 
the (100) surface of palladium (for the geometry of the 
(100) surface see the inset of Fig. ||). The employed 
high-dimensional PES was derived from density func- 
tional theory calculations Q which revealed that it is 
strongly corrugated and that many activated and only 
few non-activated paths towards dissociative adsorption 
exist. Figure m shows a cut through this PES, which 
displays one of the few non-activated pathways towards 
dissociated adsorption. In this example the II2 molecule 
approaches with an orientation parallel to the surface 
and its center of mass at the bridge position (see inset of 

Fig.g). 

The six-dimensional quantum dynamical calculations 
show that the probability for dissociative adsorption of 
H2 at the clean Pd(lOO) is very high at low kinetic 
translational energies of the impinging molecules (up 
to 70 %), and that it decreases substantially when the 
kinetic translational energy is increased (to 25 % at 
Ei =0.1 eV). This result is displayed as the full curve in 
Fig. y; it is in fact in good agreement with experimental 
findings [Q. We note that thermal energies are in the 
just noted range, i.e. typically Ei < 0.1 eV. The analysis 
of these results revealed that the decrease of the sticking 
probability with increasing kinetic translational energy is 
not due to a precursor-mediated mechanism, as was com- 
monly believed, but due to dynamical steering [^: If the 
impinging molecule is slow the attractive forces of the 
PES steer it towards a favorable (non-activated) path- 
way, even if the initial conditions were pointing towards 
an unfavorable direction. However, if the molecule is fast, 
its inertia hinders its reorientation. Such a molecule will 
then probably experience an energy barrier and thus will 
be reflected. Of course, a very fast molecule will simply 
overrun energy barriers. Another result of the calcula- 
tions was that the dissociation should be hindered if low 
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FIG. 1. Contour plot of the PES along a two-dimensional 
cut through the six-dimensional coordinate space of 
H2/Pd(100). The inset shows the orientation of the molec- 
ular axis and the lateral H2 center-of-mass coordinates. The 
coordinates in the figure are the H2 center-of-mass distance 
from the surface Z and the H-H interatomic distance d. The 
dashed line is the optimum reaction path. Energies are in eV 
per H2 molecule. The contour spacing is 0.1 eV. 



energy molecules are rotationally excited, the effect also 
being dependent on the molecular orientation, and this 
prediction has been indeed confirmed by recent experi- 
ments |ll|,0. 

In this paper we will discuss the importance of quan- 
tum effects in dissociation and adsorption. At first we 
note that calculated H2 sticking probability exhibits a 
strong oscillatory structure as a function of the inci- 
dent energy (see the dashed line in Fig. g). These os- 
cillations are a direct consequence of the quantum na- 
ture of H2 scattering. In a simplified description, i.e., 
neglecting rotational and vibrational degrees of free- 
dom of H2 , elastic scattering at a periodic surface gives 
rise to reflected beams, travelling in directions (k|| -|- 
g, —\/yi\ — 2k|| • g — g2). Here (kii,/;:^) is the wave vec- 
tor of the incident II2 beam and g is a two-dimensional 
reciprocal lattice vector of the surface. The condition for 
emerging beams is that the argument under the square 
root is positive. Just before a new beam can emerge, 
it is already built up though it remains confined to the 
surface. Thus, this beam can not (yet) be observed di- 
rectly, but, as it is coherent with the other beams, these 
others will exhibit sharp resonance structures. In ad- 
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FIG. 2. Probability for dissociative adsorption versus ki- 
netic translational energy for a H2 beam under normal in- 
cidence on a clean Pd(lOO) surface. The six-dimensional 
quantum dynamical results (solid and dashed line) are for 
molecules initially in the vibrational and rotational ground 
state. For the full line an energy spread typical for beam 
experiments has been assumed. Diamonds show results of 
a classical treatment of the H nuclei (initially non-vibrating 
molecules) and circles correspond to a quasi-classical treat- 
ment, i.e. the molecules vibrate initially with an energy of 
0.26 eV, which is the zero-point energy of H2. The inset 
shows an enlargement of the results at low-energies. 



dition, oscillations can also be caused by selective ad- 



sorption resonances |13|. These structures are known in 
He and H2 scattering since the 1930s |1J] and in elec- 
tron scattering jl^Jl^l, however, for the particular sys- 
tem H2/Pd(100) they have not been observed yet. It is 
obvious that an experimental verification requires spe- 
cial care as the strengths and energies of the oscillatory 
structures depend sensitively on the initial conditions, 
and on the quality of the surface. For example, surface 
imperfections (e.g. adatoms and steps) will destroy the 
mentioned coherence and thus suppress the oscillations. 

In previous experiments pl3| the molecular beam was 
not strictly mono-energetic. By taking the experimental 
energy spread into account, /S.E/Ei = 0.2 ||l^, the oscil- 
lations of the quantum-mechanical study are smoothed 
out and the dashed line in Fig. || is turned into the solid 
one. 

We have now complemented the six-dimensional quan- 
tum dynamical study by classical molecular dynamics 
simulations using exactly the same PES. This allows for 
the first time a thorough comparison of classical and 
quantum dynamics of H2 scattering, dissociation, and 
adsorption. 

A common approach in classical calculations aiming 
at taking some quantum effects of the H2 molecule into 
account is to give the hydrogen molecules initially a 



vibrational energy which equals the H2 zero-point en- 
ergy |6|J7|, p^ , [lq | . This is called a quasi-classical ap- 
proach. The open circles in Fig. display the corre- 
sponding results. It is quite obvious that they are quan- 
titatively and even qualitatively at variance with the cor- 
rect, i.e. quantum-dynamical, result (dashed curve in 
Fig. pi). Firstly, the quasi-classical calculations do not 
show any oscillatory structure. Secondly, they do not 
exhibit any pronounced initial decrease with increasing 
kinetic translational energy, i.e. steering is absent in the 
quasi-classical calculations. Only at very high energies, 
where quantum effects play a minor role, quantum and 
quasi-classical trajectory calculations get into agreement. 

The diamonds in Fig. g correspond to a strict classical 
treatment where the H2 molecules have initially no vibra- 
tional energy. It is interesting to note that these calcu- 
lations agree much better with the quantum-mechanical 
results than the quasi-classical one. Again, the quan- 
tum oscillations are obviously absent, but now the steer- 
ing effect is in operation. The absence of the oscillatory 
structure in the classical and quasi-classical calculations 
is caused by the fact that the coherence and interference 
of molecular beams is not present in classical mechanics. 

We will now show that further quantum effects, namely 
the generation of additional zero-point energies, lead to 
the qualitiative difference between the quasi-classical and 
the averaged quantum dynamical calculations. The in- 
teraction of the II2 molecule with the surface gives rise to 
these additional zero-point energies (due to hindered ro- 
tations and hindered surface diffusion) and this energy, 
obviously, needs to be taken from that of the incident 
beam. For the ease of the following discussion we will 
assume that the approaching molecule travels along the 
optimum reaction pathway (full line in Fig. ^). Let us 
at first describe the correct, i.e. quantum-mechanical, 
results. When the molecule gets close to the surface, 
the energy of its stretch vibrations is lowered, because 
when the H-surface bond starts to build up, the H-H 
bond strength will at the same time soften. Since the vi- 
brational frequency is high, the change of the frequency 
during one oscillation is small compared to the frequency, 
i.e. 
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This suggests that the H2 vibration follows the molec- 
ular center of mass motion adiabatically, which is in- 
deed confirmed within a very good approximation by 
the quantum mechanical M] as well as by the classi- 
cal calculations. Interestingly, the change in zero-point 
energy is not transferred to translational motion. Thus, 
in contrast to common believe, it cannot be utilized to 
overcome energy barriers. The reason is that at the 
point where the H-H stretch-frequency is lowered, i.e. 
where the H-H bond is already weakened, also a notice- 
able H-surface bond has already started to develop, as 



mode 


zero-point energy (eV) 


H-H vibration 

polar rotation 

azimuthal rotation 

translation perpendicular to bridge 

translation parallel to bridge 


0.126 
0.084 
0.033 
0.051 
0.056 


sum 


0.350 



TABLE I. Zero-point energies of the H2 molecule at the 
point where the sum of all zero-point energies has its maxi- 
mum (marked by X in Fig. hi) . The H2 center-of mass is above 
the bridge position between two adjacent Pd atoms (see inset 
of Fig. hi). The energies are given per H2 molecule. 



was noted already above. In other words, close to the 
surface the PES is strongly corrugated. This leads to 
the building up of zero-point energies in the degrees of 
freedom of the hydrogen molecule perpendicular to the 
optimum reaction path. In table | all these calculated 
zero-point energies are listed at the point where their 
sum has its maximum. This geometry is indicated in 
Fig. n^ by the cross. Although the zero-point energy of 
the stretch vibrations is significantly lowered, this lower- 
ing is over-compensated by the zero-point energies in the 
other modes so that the sum of all zero-point energies, 
0.35 eV per H-H pair is even larger than the H2 gas- 
phase value of 0.26 eV. That the sum of all zero-points 
energies in the case of hydrogen dissociative adsorption 
is approximately equal or even higher than in the gas- 
phase has also been found in total-energy calculations 
for H2/A1(110) |§ and H2/Si(100) |l|, where values of 
0.25 eV/molecule and 0.29 eV/molecule have been deter- 
mined at the top of the barrier to dissociative adsorption. 

Obviously, in classical or quasi-classical calculations 
the concept of zero-point vibrations does not exist. 
Therefore in the quasi-classical calculations the vibra- 
tional energy can be almost fully transferred to the trans- 
lational energy which enables the molecule to traverse 
energy barriers with a height larger than its initial ki- 
netic translational energy. This is actually the reason 
why the steering effect is absent in the quasi-classical re- 
sults (open dots in Fig. 0): These vibrating molecules 
do not experience a strong energy barrier because the vi- 
brational energy is efficiently transferred to translational 
energy. 

The analysis of the zero-point energies along the opti- 
mum reaction pathway (see Fig. |^) teaches us that the 
sum of all zero-point energies is roughly constant during 
the dissociation. This approximate "constant of motion" 
can be simulated in classical trajectory calculations by 
ignoring the zero-point energies all along the way. As 
mentioned above, this strictly classical treatment (see the 
diamonds in Fig. g) gives results in remarkable agreement 



with the averaged quantum dynamical results at low and 
high energies. The disagreement in the medium energy 
regime could perhaps be due to the fact that the sum of 
the zero-point energies is not really constant but becomes 
larger than the gas-phase zero-point energy so that the 
quantum mechanical sticking probability is lower due to 
the effectively higher barriers. 

Our results clearly show that for kinetic translational 
energies between 0.15 and 0.6 eV the probability for dis- 
sociative adsorption is reduced by quantum effects by 
50 % compared to the quasi-classical result. Thus, in 
contrast to the common belief, quantum effects in the 
dynamics of H2 chemical reactions at surfaces are clearly 
not negligible. 

An influence of the zero-point energies in traversing 
barrier regions has also been found in the comparison of 
three-dimensional quantum and classical calculations of 
the dissociative adsorption of H2/Cu(lll) |^. However, 
since these calculations were restricted to low dimensions, 
the difference between quantum and classical results was 
much less pronounced because the sum of the zero-point 
energies was smaller. From the above discussion it is 
clear that a proper theoretical description of the role of 
zero-point effects requires a high-dimensional quantum- 
mechanical treatment of the nuclear motion. 

The problem of a proper treatment of zero-point en- 
ergies in quasiclassical trajectory calculations is well- 
known, especially in the gas-phase community |17| , p^ . 
One possible way to deal with this problem is the re- 
duced dimensionality treatment in the vibrationally adi- 
abatic approximation (for a overview see Ref. p2[). In 
this approach a small number of degrees of freedom is 
treated dynamically while the remaining degrees of free- 
dom are taken into account by adding the sum of their 
zero-point energies to the potential along the reaction 
path. Another more elaborate approach is to constrain 
the energy in each vibrational mode to be greater than 
its zero-point value [ [l7|Jlq| . 

In our purely classical approach we ignore zero-point 
energies all along the reaction path. But this approach 
is actually in the spirit of the vibrationally adiabatic ap- 
proximation. It effectively takes the zero-point energies 
into account through a shift of the potential along the 
reaction path corresponding to the sum of all zero-point 
energies. This shift, however, is constant along the reac- 
tion path. Moreover, we still keep the full dimensionality 
of the problem by explicitly treating all degrees of free- 
dom dynamically. This is indeed essential since for ex- 
ample the steering effect is absent in a low-dimensional 
treatment of the Il2/Pd(100) system [E3[ (see also the 
discussion below). 

In order to illustrate some of the dynamical effects dis- 
cussed above we have plotted in Figs. and H some typi- 
cal trajectories of the classical molecular dynamics calcu- 
lations. The upper panel shows snapshots of trajectories 
of the impinging II2 molecules while the lower panel dis- 
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FIG. 3. Classical trajectories of a hydrogen molecule im- 
pinging on a Pd(lOO) surface. The initial conditions are 
chosen in such a way that the trajectories are restricted 
to the xz-plane. Upper panel: snapshots of the trajecto- 
ries; lower panel: trajectories in the xz-plane. Trajectory 
A: quasi-classical trajectory with an initial kinetic transla- 
tional energy of 0.12 eV. Trajectory B: same initial condi- 
tions as in trajectory A, except that the molecule is initially 
non-vibrating. Trajectory C: same initial conditions as in 
trajectory B, except that the molecule has a lower kinetic 
translational energy of 0.01 eV. 



plays the whole trajectories. The initial conditions are 
chosen in such a way that the trajectories are restricted 
to a xz-plsme. Trajectories A and B in Fig. illustrate 
the difference between quasi-classical and classical cal- 
culations. They have identical initial conditions with a 
kinetic translational energy of 0.12 eV except for the fact 
that in case A the molecule initially vibrates with an en- 
ergy corresponding to the zero-point energy while in case 
B the molecule is initially non- vibrating (due to the small 
amplitude of 0.12 A the vibrations are hardly perceivable 
on the scale of the figure). For the chosen initial condi- 
tions of the impinging molecules there is an energy bar- 
rier hindereing the dissociative adsorption, and the initial 
translational energy is not large enough for the molecule 
to overcome this. Thus, the molecule is scattered back 
at the repulsive potential (trajectory B), whereas if the 
initial vibrational zero-point energy is transferred to the 
translation due to the adiabaticity of the vibrations (tra- 
jectory A), the molecule can directly dissociate. 

Finally, trajectory C illustrates the steering effect [g,||. 
The initial conditions are as in trajectory B, but now the 
molecule has a much lower kinetic translational energy, 
Ei = 0.01 eV. The molecule is so slow that the attractive 
forces can reorient the molecule so that it can follow a 
non-activated path towards dissociative adsorption. We 
like to emphasize that steering works in all geometry pa- 
rameters of the nuclei and becomes more important when 



the configuration space is high-dimensional - in a one- 
dimensional description (just restricting the reaction to 
"the optimum" reaction pathway) or a two-dimensional 
description (restricting the dynamics to only one elbow 
p3[), the steering effect is absent. 

Trajectory C in Fig. ^ represents a classical illustration 
of steering. We note that this picture also applies when 
the quantum nature of H2 molecules is considered, and 
we now briefly address the quantum mechanical picture. 
Let us assume that the incident H2 molecule is in the 
rotational ground state, ji = (the description for ro- 
tationally excited states is analogous). The ji — state 
refers to an isotropic distribution of molecular orienta- 
tion. Thus, the molecule is not rotating, but when the 
orientations of different molecules in the beam were mea- 
sured, one would find all orientations with the same prob- 
ability. In other words, the probability that a molecule 
impinges with the favorable orientation noted in Fig. 1 
is very low. For molecules with low kinetic energy the 
steering effect will align all molecules, and trajectory C 
shows one example, which acts on the molecular polar 
coordinate 6. Equivalent processes take place in the az- 
imuthal cj) and the {X, Y) coordinates. Also with re- 
spect to the {X, Y) coordinates of the H2 molecule we 
note their meaning in (basic) quantum mechanics: A 
plane wave of hydrogen molecules representing a molecu- 
lar beam impinging under normal incidence on a surface 
has no {X, Y) dependence in the gas-phase. This means 
that the probability to find a molecule at any given point 
point in a {X, Y) plane far away from the surface is the 
same. This brief reminder to the concepts of quantum 
mechanics highlights also the big advantage of quantum 
dynamics in comparison to molecular dynamics. A quan- 
tum dynamical approach includes the ensemble of all H2 
orientations {9, cf)) and all (X, Y) points of impact auto- 
matically. In a molecular dynamics calculations the same 
ensemble has to be built up step by step, i.e. by perform- 
ing thousands of molecular dynamics runs, namely one 
run for each choice of {X, Y, 9, (j)). 

The trajectory in Fig. ^ finally indicates a process 
which in a quantum-mechanical system will give rise to 
an oscillatory structure (see Fig. g). The initial condi- 
tions are as in the trajectory of Fig. pC, only the impact 
point in the surface unit cell and the inital polar orienta- 
tion are slightly changed. The snapshots of the trajectory 
in the upper panel are taken approximately every 150 fs. 
They are consecutively numbered so that the dynamics 
can be roughly followed. Due to the strong anisotropy 
and the corrugation of the potential the translational en- 
ergy perpendicular to the surface is efficiently transferred 
to rotational and kinetic energy parallel to the surface. 
Hence the molecule is temporarily trapped; in this case 
the molecule spends more than 5 ps in front of the sur- 
face before it is dissociated. Note that this trapping does 
not occur because of energy transfer to the surface but 
because of the energy redistribution among the II2 de- 
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FIG. 4. Classical trajectory of a hydrogen molecule imping- 
ing on a Pd(lOO) surface. The initial kinetic energy is 0.01 eV, 
the initial conditions are only slightly different from those of 
Fig. He (see text). Upper panel: snapshots of the trajectories 
which are taken approximately every 150 fs and consecutively 
numbered. Lower panel: trajectories in the xz-plane; for the 
sake of clarity the trajectories are shifted in the 2-direction 
with respect to each other. 



grees of freedom. Quantum mechanically the trapped 
molecules correspond to metastable states at the surface 
which only exist for certain discrete energies; in addition, 
only discrete values for the energy and momentum trans- 
fer are allowed. Thus only for particular initial energies 
the dynamical trapping into these so-called selective ad- 
sorption resonances is possible [O 14 1. As mentioned in 
the discussion of Fig. 0, these resonances in addition to 
the opening up of scattering channels lead to the sharp 
peaks in the quantum-mechanical sticking probability. 

In conclusion, we have shown that there are large dif- 
ferences between classical and quantum dynamics of hy- 
drogen. The classical results do not show any oscilla- 
tory structure which in the quantum dynamics is due 
to the opening up of new scattering channels and res- 
onances. Even more remarkable, zero-point effects can 
cause substantial deviations between classical and aver- 
aged quantum dynamical calculations. A corrugated and 
anisotropic potential energy surface leads to the building 
up of zero-point energies which effectively increase the 
minimum potential in the quantum calculations. This 
can change the dynamics of chemical reactions not only 
quantitatively, but also qualitatively, and thus strongly 
affects the rate constants of catalytic reactions involving 
hydrogen. These zero-point effects will also be relevant 
for other studies of hydrogen dynamics in corrugated po- 
tentials, like for example in hydrogen diffusion on surfaces 
or in the bulk. If, however, the sum of the zero-point en- 



ergies is approximately constant during the traversing of [23] G. R. Darling and S. HoUoway, Surf. Sci. 268, L305 
regions with energy barriers, which for hydrogen disso- (1992). 

elation on surfaces seems to be the case for a wide class 
of systems, the discrepancies between classical and aver- 
aged quantum dynamics can be diminished by neglecting 
zero-point energies in the initial conditions. 
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